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Abstract.  Rarefied  gas  flows  in  ultra- thin  film  slider  bearings  are  studied  in  a  wide  range  of  Knudsen  numbers.  The 
generalized  Reynolds  equation  first  derived  by  Fukui  and  Kaneko  [1],  [2]  has  been  extended  by  allowing  for  bounding  surfaces 
with  different  physical  structures,  as  an  issue  of  relevance  for  applications.  Since  the  solution  of  this  equation  requires  that  the 
Poiseuille  and  Couette  flow  rates  between  two  parallel  plates  have  to  be  accurately  calculated  in  advance,  we  have  used  our 
recent  results  on  Poiseuille  flow  [3],  [4]  and  new  results  on  Couette  flow  to  evaluate  the  lubrication  characteristics. 


FLOW  OF  A  RAREFIED  GAS  BETWEEN  TWO  PARALLEL  PLATES:  THE 
POISEUILLE-COUETTE  PROBLEM 

Let  us  consider  two  plates  separated  by  a  distance  h  and  a  gas  flowing  parallel  to  them,  in  the  x  direction,  due  to 
a  pressure  gradient.  The  lower  boundary  (placed  at  z  =  —h/2)  moves  to  the  right  with  velocity  U ,  while  the  upper 
boundary  (placed  at  z  =  h/2)  is  fixed.  Both  boundaries  are  held  at  a  constant  temperature  T0. 

If  the  pressure  gradient  is  taken  to  be  small  as  well  as  the  velocity  U,  it  can  be  assumed  that  the  velocity  distribution 
of  the  flow  is  nearly  the  same  as  that  occurring  in  an  equilibrium  state.  This  means  that  the  Boltzmann  equation  can 
be  linearized  about  a  Maxwellian  /0  by  putting  [5]: 


f  =  f0(l+h)  (1) 

where  /(jt,z,c)  is  the  distribution  function  for  the  molecular  velocity  c  expressed  in  units  of  (2 ,RT0)1/2  ( R  being  the 
gas  constant),  z  is  the  coordinate  normal  to  the  plates  and  h(z,  c)  is  the  small  perturbation  upon  the  basic  equilibrium 
state.  If  one  assumes  the  linearized  BGK  model  for  the  collision  operator  [6],  the  Boltzmann  equation  reads  [7]: 


where  by  definition 
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with  p  and  p  being  the  gas  pressure  and  density,  respectively,  and  6  is  the  collision  time.  Consequently,  the  integral 
equation  for  the  bulk  velocity  of  the  gas  can  be  written  as  follows: 

q(z)  =  7l~?  t  e~°zi Z(z,cZl)dcZi  (3) 

From  Eq.  (2)  we  obtain  in  integral  form  [3],  [8]: 
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Z(z,cz)  =  exp(-(z+|sgncz)/(cz0))Z(-|sgncz,cz)  + 

[  h  [q(t)-kO/2]/(czO)dt 

J~jsgncz  \  \CZ\V  J 

with  the  values  at  the  boundary,  Z(— ^sgncz,cz),  depending  on  the  model  of  boundary  conditions  chosen.  In  the 
following,  we  will  consider  the  Maxwell  boundary  conditions  and  specialize  the  analysis  to  walls  having  different 
physical  properties  so  that  two  accommodation  coefficients  ( al:a2 )  can  be  defined.  In  this  case,  the  boundary 
conditions  can  be  written  as  [3],  [8]: 


Z+(h/2,cz)  =  (1  -ax)Z-(h/2,-cz)  (5) 

Z+(-h/2,cz)  =  a2U  +  (l-oc2)Z-(-h/2,-cz)  (6) 

where  U  is  expressed  in  units  of  (2RT0)]I2\  Z~(—h/2,cz),  Z~{h/2,cz)  are  the  distribution  functions  of  the  molecules 
impinging  upon  the  walls  and  Z+(— /i/2,cz),  Z+(/z/2,cz)  the  distribution  functions  of  the  molecules  reemerging  from 
them. 

Once  the  function  at  the  boundary,  Z(— ^sgncz,cz),  has  been  evaluated  following  the  analytical  procedure  reported 
in  [3],  [4],  the  substitution  of  the  integral  formula  (4)  in  the  definition  (3)  of  q(z)  gives  the  integral  equation  for  the 
bulk  velocity  of  the  gas: 


q(z)  =  \kS[  1  -  \l/p(u)]  +  U yc(u)  (7) 

Eq.  (7)  shows  that  the  gas  velocity  is  induced  by  the  superposition  of  two  distinct  effects.  The  gas  moves  by  an  imposed 
pressure  gradient  (Poiseuille  flow)  and  by  the  shear  driven  flow  due  to  the  motion  of  the  bottom  surface  (Couette  flow). 
The  non-dimensional  functions  \\fp{u)  and  i j/c(u),  giving  the  Poiseuille  and  Couette  contributions,  respectively,  can  be 
written  as: 


i  rs/ 2  { 
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where  Tn  (v)  is  the  Abramowitz  function  defined  by 

r+™ 

Tn(x)  =  /  tnexp(—t2—x/t)dt 

Jo 

Sn(x)  is  a  generalized  Abramowitz  function  defined  by 

r+oo  f”exp(— t2  —  x/t) 


n  (-^j  ^  5  5  ^2  )  y 


■  (1  —  oq)(l  —  a2)exp(— 2  8/t) 
and  the  following  non-dimensional  variables  have  been  introduced: 


dt 


(9) 


5  =  h/ 6.  w  =  t/07  u  =  z/0. 
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FIGURE  1.  Geometry  of  a  slider  bearing. 


Using  Eq.  (7),  the  flow  rate  (per  unit  time  through  unit  thickness)  defined  by  [7]: 


can  be  expressed  as  the  sum  of  the  Poiseuille  flow  (Fp)  and  the  Couette  flow  (Fc)  as  follows: 

F  =  Fp  +  Fc  =  —  Qp{&i  7  ^2)  3  2  2c (5,  oq ,  oq) 

where 


1  1  r5/2 

Qp(S,ccvcc2)  =  --  +  j2  j  \ifp(u)du 

2  f8/2 

Qc(8,ava2)  =  -  /  yc(u)du 
o  J-8/2 

are  the  non-dimensional  volume  flow  rates. 


(10) 


(11) 


THE  GENERALIZED  REYNOLDS  EQUATION 


The  analysis  developed  in  the  previous  section  can  be  applied  to  the  slider  bearing  problem  in  lubrication  theory. 
The  basic  geometry  of  the  two-dimensional  gas  film  is  outlined  in  Fig.  1.  Unlike  the  configuration  considered  in  the 
previous  section,  the  upper  plate  is  slightly  inclined  at  a  small  angle  7.  Since  the  pitch  angle  7  is  typically  less  than 
1°,  the  pressure  p  in  the  gas  is  taken  to  be  constant  in  the  vertical  direction,  while  at  the  left  (. x  =  0)  and  right  (x  =  L) 
boundaries  it  is  fixed  at  ambient  pressure  pQ . 

In  standard  hydrodynamic  lubrication  theory,  the  pressure  distribution  in  the  gas  film  can  be  computed  from  the 
classical  Reynolds  equation  (an  approximation  of  the  full  Navier-Stokes  equation): 


—  (pH3—-APh)  =0  (12) 

dX\dXj  K  J 

where  the  no-slip  boundary  conditions  for  the  velocities  at  the  bottom  and  top  surfaces  have  been  used.  Equation  (12) 
is  written  in  terms  of  the  following  dimensionless  quantities 


and  the  bearing  number  is  defined  as 


P  = 


H  = 


Po 


(13) 


_  6 nUl 

Poh2a 


(14) 


where  p  is  the  dynamic  viscosity  of  the  gas. 

To  extend  the  generalized  Reynolds  equation,  the  mass  flow  conservation  equation  across  the  film  thickness  must 
be  applied  [1],  [2],  with  the  flow  rate  F  given  by  Eq.  (11): 
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(15) 


Tx^hlQp^iav°^~^Qc^,ava^)  =° 

For  the  purpose  of  a  direct  comparison  with  the  classical  Reynolds  equation  (12),  let  us  introduce  the  Poiseuille  relative 
flow  rate  [1],  [2]: 


Qp(S,ava2) 


Qpi^a^c^) 

Qcon 


(16) 


where  Qcon  =  8/6.  Furthermore,  the  rarefaction  parameter  5  can  be  expressed  as:  8  =  80PH ,  where  80  is  the 
characteristic  inverse  Knudsen  number  defined  by  the  minimum  film  thickness,  h0,  and  the  ambient  pressure  p0  as: 


s  _  PoK 

°~  nV^RTo 

Finally,  assuming  that  the  heat  generation  in  the  gas  is  very  small,  so  that  an  isothermal  process  can  be  considered,  the 
non-dimensional  generalized  Reynolds  equation  reads: 

^  (Qp(80PH,ava2)PH3^-Qc(S0PH,al,a2)APHSj  =0  (17) 

If  the  two  walls  are  identical  (oq  =  a2  =  a),  the  Couette  flow  rate  is  independent  of  the  Knudsen  number  regardless 
of  the  value  of  the  accommodation  coefficient  a  : 


Qc(80PH,a)  =  1  Va,  V<5. 

Therefore,  in  this  case,  Eq.  (17)  reduces  to  the  generalized  Reynolds  equation  introduced  by  Fukui  and  Kaneko  [1], 
[2],  while  in  the  continuum  regime,  Eq.  (17)  reduces  to  the  classical  Reynolds  equation  (12). 

Writing  the  non-dimensional  film  thickness  H  in  terms  of  the  longitudinal  coordinate  X , 


such  that 


/ 

L 


dP 

dX 


Eq.  (17)  can  be  analytically  integrated  to  give: 


IP 

L  \hB  )  dH 


(18) 
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,3  dP 


Qp(80PH,  a, ,  a^PH3  —  +  QC{80PH,  avoc2)APH  =  Kx 

drl 


where  Kx  is  a  constant  of  integration.  The  substitution  of 


(19) 


in  Eq.  (19)  gives: 


PH  =  t; 


(20) 


_  _C _ [Qc(80^  oq,  oc2)A^  K/\ 

dH  ~  H  l/L(h2/h0-l)Qp(80^ava2)Hi ; 

Eq.  (21)  has  been  solved  numerically  using  the  relaxation  methods.  To  apply  this  numerical  scheme,  the  differential 
equations  have  to  be  replaced  by  finite-difference  equations  on  a  point  mesh.  The  solution  of  the  resulting  set  of 
equations  is  determined  by  starting  with  a  guess  and  improving  it  iteratively  using  Newton’s  method.  The  Poiseuille 
flow  rate  coefficient  Qp(8,al,a2)  has  been  evaluated  by  means  of  the  numerical  method  described  in  [3]  and  a 
variational  technique  which  applies  to  the  integrodifferential  form  of  the  Boltzmann  equation  based  on  the  BGK 
model  [4].  In  order  to  compute  the  Couette  flow  rate  Qc(S,a1:a2)  we  have  solved  numerically  Eq.  (9),  extending  a 
finite  difference  technique  first  introduced  by  Cercignani  and  Daneri  [7] . 
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FIGURE  2.  Pressure  profile  for  80  —  0.5.  The  line  styles  indicate  a  —  0.8  (solid),  a  —  0.3  (dashed),  and  a  —  0.1  (dot  dashed). 
From  left  to  right,  the  bearing  number  A  is  10,  50,  200. 


FIGURE  3.  Pressure  profile  for  80  —  0.5.  The  line  styles  indicate  a{  —  0.5  a2  =  0.8  (solid),  ax  —  0.5  a2  =  0.3  (dashed),  and 
ax  =  0.5  a2  =  0. 1  (dot  dashed).  From  left  to  right,  the  bearing  number  A  is  10,  50,  200. 


NUMERICAL  RESULTS  OF  LUBRICATION  CHARACTERISTICS 

Once  £(//)  has  been  numerically  evaluated  on  a  grid  that  spans  the  domain  of  interest,  Eqs.  (18)  and  (20)  give  the 
pressure  field  in  the  gas  film  as  a  function  of  X.  Furthermore,  a  prediction  of  the  vertical  force  acting  on  the  upper 
surface  of  the  slider  bearing,  crucial  for  practical  design,  may  be  obtained  from  the  load  carrying  capacity  W ,  defined 
as 


/  rL/1 

W=-  (P-l)dX  (22) 

L  Jo 

In  order  to  investigate  the  effects  of  the  rarefaction  parameter  80  and  the  bearing  number  A  on  the  basic  lubrication 
characteristics  (pressure  distribution  and  load  carrying  capacity),  the  parameters  describing  the  gas  film  geometric 
configuration  were  fixed  at  the  following  values:  h2/h0  =  2,  L/h0  —  100.  Figure  2  shows  the  pressure  field  as  a 
function  of  the  longitudinal  coordinate  X  at  three  different  bearing  numbers:  A  =  10,  50,  200.  To  assess  the  influence 
of  the  boundary  conditions,  the  profiles  corresponding  to  different  accommodation  coefficients  (for  bounding  walls 
supposed  physically  identical)  are  drawn  in  each  panel  of  Fig.  2.  Looking  at  the  pictures,  one  sees  that  the  pressure 
distribution  in  the  gas  film  increases  with  increasing  A.  Furthermore,  at  fixed  bearing  number,  the  pressure  field  reduces 
by  increasing  the  fraction  of  gas  molecules  specularly  reflected  by  the  walls.  Fig.  3  reports  the  pressure  profiles  for 
the  same  parameters  as  in  Fig.  2,  except  that  now  the  two  bounding  plates  are  allowed  to  re-emit  the  impinging  gas 
molecules  differently,  so  that  two  accommodation  coefficients  can  be  defined.  We  keep  the  accommodation  coefficient 
of  the  upper  wall  (cq)  fixed  and  vary  the  other  one  (a2).  A  comparison  with  Fig. 2  shows  that,  for  every  A,  the  pressure 
distribution  significantly  depends  on  a2  and  only  weakly  on  .  This  picture  remains  true  at  different  Knudsen  numbers 
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FIGURE  4.  Pressure  profile  for  A  =  50.  The  line  styles  indicate  a  —  0.8  (dot  dashed),  ax  —  0.5  a2  —  0.8  (long  dashed),  a  —  0.3 
(solid),  and  ax  —  0.5  a2  —  0.3  (dashed).  The  inverse  Knudsen  number  80  is  10”3  (left)  and  10  (right). 


FIGURE  5.  Pressure  profile  versus  X.  Comparison  between  our  results  (solid  line)  and  DSMC  data  [10]  (open  circles).  From  left 
to  right,  the  parameters  are:  (a)  80  =  0.7,  A  =  61.6,  a  =  1;  (b)  S0  =  0.2,  A  =  1264,  a  —  1;  (c)  80  =  0.7,  A  =  61.6,  a  =  0.7. 


progressing  from  free  molecular,  through  transitional,  to  continuum  regions  (see  Fig.4).  It  is  worth  noting  that,  when  A 
increases,  the  Couette  contribution  to  the  lubrication  flow  rate  becomes  dominant  compared  with  the  Poiseuille  flow. 
Therefore,  if  the  two  walls  are  identical,  the  influence  of  the  Knudsen  number  on  the  load  carrying  capacity  decreases 
as  A  increases,  since  Qc  is  independent  of  8  and  a.  On  the  contrary,  if  the  two  walls  have  a  different  physical  structure 
the  load  carrying  capacity  shows  a  dependence  on  both  the  Knudsen  number  and  the  accommodation  coefficients  ax , 
a2.  For  the  validation  of  the  code,  our  outputs  have  been  compared  with  the  results  from  DSMC  (Direct  Simulation 
Monte  Carlo)  simulations  published  by  Alexander  et  al.  [10]  in  the  case  of  Maxwell’s  boundary  conditions  on  two 
physically  identical  walls  (see  Fig.5). 
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